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O ■ ABSTRACT 

(N 

S ■ Using the Spitzer Space Telescope and the VLA, we present the first direct evi- 

dence that the well-known far-infrared/radio correlation is valid to cosmologically 
significant redshift. We also confirm, with improved statistics compared with pre- 
vious surveys, a similar result for the Mid-IR/radio correlation. We explore the 
dependence of monochromatic ^24 and qm on z. The results were obtained by 
cn . matching Spitzer sources at 24 and 70//m with VLA 1.4GHz /iJy radio sources 

v^ ■ obtained for the Spitzer FLS. Spectroscopic redshifts have been obtained for over 

Q . 500 matched IR/radio sources using observations at WIYN, Keck and archival 

O ■ SDSS data extending out to z > 2. We find that q24 shows significantly more 

dispersion than q^q. By comparing the observed fiuxes at 70, 24 and 4.5/xm with 
a library of SED templates, we find that the larger dispersion in 524 is predictable 
in terms of systematic variations in SED shape throughout the population. Al- 
To . though the models are not able to encompass the full range of observed behavior 

(both the presence of either extremely fiat or extremely steep IR SEDs), the fit- 
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ting parameters were used to "k-correct" the higher-z galaxies which resulted in a 
reduced scatter in q. For comparison, we also corrected these data using the SED 
for M82. The results for 24 and 70/im provide strong consistent evidence for the 
universality of the mid-IR/radio and far-IR/radio correlations out to redshifts of 
at least z = 1. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: 
high-redshift — galaxies: IR/radio correlation 
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1. Introduction 

It has now been more than a quarter of a century since the first hints emerged of a very 
tight correlation between the global far-infrared (FIR) and radio emission from galaxies. 
Initially based on comparisons between challenging ground-based lO/xm/radio correlations 
of small samples of galaxies (van der Kruit 1973; Condon et al. 1982; Rikard & Harvey 1984), 
the remarkable nature of the FIR/radio correlation became increasingly apparent with larger 
IR-samples obtained from the Infrared Astronomical Satellite (IRAS) all sky survey (Dickey 
& Salpeter 1984 ; De Jong et al. 1985; Helou, Soifer, & Rowan-Robinson 1985; Condon 
& Broderick 1986; Wunderlich et al. 1987; Hummel et al. 1988; Fitt, Alexander, & Cox 
1988). IRAS showed that the FIR-radio correlation appeared to apply to a wide range of 
Hubble-types, extending over more than three orders of magnitude in IR luminosity from 
dust-rich dwarfs to ultra-luminous infrared galaxies (ULIRGS). 

Though revolutionary, IRAS could only systematically study the relatively nearby uni- 
verse. The Infrared Space Observatory (ISO) probed, for the first time, galaxy populations 
in the higher-z universe, but had limited long-wavelength coverage and sensitivity. Based 
on the deepest ISO fields, it was evident that even at 15;um there was a loose correlation 
between the Mid-IR (MIR) emission, and the radio continuum (Cohen et al. 2000, Elbaz 
et al. 2002, Gruppioni et al. 2003). Using M82 as a template, Garrett (2002) extrapolated 
from 15;um to longer wavelengths, and suggested that the Far-IR/radio correlation proba- 
bly extends to 2; > 1. At the other extreme, sub- mm observations using instruments like 
SCUBA, began to reveal a new population of likely star-formation active galaxies. Although 
the statistics and redshift determinations are still growing (see Chapman et al. 2003), the 
correlation between /iJy radio and SCUBA sources suggest that the FIR/radio correlation 
is likely to be valid at much higher redshift (see Carilli & Yun 1999, Ivison et al. 1998, 
Frayer et al. 1999, Ivinson et al. 2002). The universality of the correlation, if confirmed, is 
important because it implies that galaxies observed at large look-back times share many of 
the same properties as fully-formed mature galaxies seen locally. 



2. Spitzer, VLA Radio and Spectroscopy Observations 

Spitzer observations^ of the First Look Survey (FLS) were obtained from 1-11 December 
2003 (Storrie-Lombardi et al. in prep., Faddaet al. 2004). A total of 4 sq. degree was covered 



^This work is based [in part] on observations made with the Spitzer Space Telescope, which is operated 
by the Jet Propulsion Laboratory, California Institute of Technology under NASA contract 1407. 
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by both the IRAC (3.6, 4.5, 5.8, 8.0/im; Fazio et al. 2004-this volume) and MIPS (24, 70 
and 160/im; Rieke et al. 2004-this volume) instruments at medium depth (5 x 12s IRAC, 
2 X 40s for MIPS 24, 1 x40s 70^m, and 2x4s at 160/im), and a smaller "verification strip" 
(0.25 square degree) was observed to an even greater depth (4 x 40s). The IRAC and MIPS 
24/im data were processed through standard SSC pipelines (Lacy, Fadda & Frayer, in prep.). 
The average FWHM of the point-spread-function of the final maps were ~2" for the IRAC 
bands, and 6, 18 and 40" for MIPS 24, 70/im and 160/xm respectively. 

IRAC sources were extracted (Lacy et al. in preparation) using the package SExtractor 
(Bertin & Arnouts 1996). MIPS 24 and MIPS 70/im sources were extracted from the images 
using the StarFinder method (Diolaiti et al. 2000). The formal 5-a sensitivities are measured 
to be 6, 7, 45 and 32/iJy for IRAC (3.6, 4.5, 5.8 & 8/im) and 0.3, 30 and 100 mJy for MIPS 
(24,70 and 160/im respectively). Tests performed on data from both the full survey and 
the verification strip at 24/im and 70/im indicate that the source flux densities show good 
reliability and repeatability down to a flux density of ~0.5 mJy and 30 mJy respectively. 

Radio observations were carried out by Condon et al.(2003) down to a formal 5-cr depth 
of 115 mJy with 35-overlapping B-array VLA pointings (FWHM of restored synthesized 
beam = 5"), revealing over 3500 /iJy radio sources. This conservative limit is appropriate 
for the entire survey, including the edges of the field which were not surveyed as deeply as 
the inner regions. Given the coverage of the MIPS and IRAC fields, which did not extend 
to the outer edges of the VLA survey, we used a slightly deeper catalog (S(20cm)> 90/iJy 
prepared by JJC) more appropriate for cross- matching with the Spitzer data. The deeper 
VLA catalog was found to provide a high degree of reliability when cross-correlated with 
deep R-band imaging of the same field. 

Our sample consists of VLA radio sources with MIPS 24 and/or MIPS 70 counterparts 
and known spectroscopic redshifts. Source matching between the Spitzer and radio catalogs 
involved an automated catalog search, followed by a visual inspection of the IR and radio 
images (this was important because of the possible presence of noise artifacts in the 70/im im- 
ages near the edges of scan legs). For the IRAC and MIPS 24/im source catalogs, we searched 
for IR counterparts to VLA sources within a radius of r < 2.5" of the VLA centroid. For the 
MIPS 70/im images, the search was performed out to a radial distance of 12" because of the 
larger effective beam-size of the 70/im observations. For MIPS 24/im and 70/im , 508 and 
227, matched sources were found respectively. Of the > 3000 matched IRAC/VLA sources, 
412 4.5/im sources were in common with the MIPS/VLA sample used in this paper^. So few 



^The reason we used IRAC 4.5/ini sources in this paper is because we have strong evidence in our sample 
that the 5.8 and 8.0/im sources are strongly evolving in flux in the z range to 0.8: presumably due to 
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160;Um detections were made that we restricted ourselves to the shorter wavelengths. 

Spectroscopy of likely MIPS sources was performed before the launch of Spitzer by 
targeting VLA radio sources in two separate surveys. The first comprised of ~ 1200 ra- 
dio sources, and were observed using the NOAO-WIYN/Hydra fiber-system (Marleau et 
al. in preparation) while the second survey yielded 80 radio sources targeted with the 
Keckll/Deimos spectrograph. The latter was part of a larger spectroscopy study of the FLS 
(Choi et al. in preparation). These data were supplemented with archival redshifts from the 
Sloan Digital Sky Survey. 



3. Model-dependent "/^-corrections" 

In this paper we present primarily flux-density ratios between the IR and the radio 
observations. Neglecting the (1 + 2;)^^ bandwidth compression term of the frequency units 
of flux density (Weedman 1986, Hogg et al. 2002), which cancels for flux ratios, the radio 
flux densities were boosted by a factor (1 + 2;)°'^, where we assumed a synchrotron power law 
of the form S oc z/~°-^, a spectral index = -1-0.7 (a value typical of an average steep spectrum 
radio source). The boosting corrects the observed flux density at the emitted frequency to 
the value it would have at the observed frequency in its rest frame. 

We approach the IR /c-correction using two methods. Where 4.5, 24 and 70/im data 
were available, we compared these data with a set of redshifted spectral energy distributions 
(SEDs) from Dale & Helou (2001; hereafter DH) and used these to compute the correction. 
As we shall see, the models only partly encompass the observed galaxy behavior. There- 
fore, as a second approach (and to extend the SED corrections to cover a larger number 
of 24/im sources), we also use a complete SED for M82 (Fadda et al. 2002) appropriately 
convolved to the resolution of the Spitzer filters to correct these data. 

The DH SEDs are a set of galaxy SEDs ranging from optical to sub-mm wavelengths 
and covering 64 different levels of excitation, as measured by the dust excitation parameter 
a ranging from 0.063 to 4.^. As discussed in detail by DH, values of a less than 1 correspond 
to a luminous starburst, while larger a values are more appropriate for describing quiescent 



the shifting out of these bands of the 3.6, 6.2 and 7.7^m(PAH) infrared features. Using the 4.5^m sources 
also increased the number of MIPS/VLA matches compared with the less sensitive IRAC filters at longer 
wavelengths 

■^The galaxy SEDare created by combining individual' "model" SEDs according to a power-law distribu- 
tion in UdM{U) ex U^^dU, where M{U) represents the dust mass heated by radiation of energy density 
U. 
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Fig. 1. — a) The 20cm radio and 70/im IR luminosity correlation (see text) for the FLS 
long-wavelength sample, b) The distribution of monochromatic Qiq values with redshift, 
uncorrected and, c) k-corrected using the DH SED-fitting method described in the text for 
the IR flux densities and assuming a (l+z)°'^ k-corrected (boosting) of the 1.4 GHz values. 

systems. To select the best suited template SED, we computed the 4.5, 24 and 70yum fluxes of 
each template redshifted to the value appropriate for each galaxy. We then use the observed 
4.5, 24 and 70yum flux densities (or 4.5 and 24|Um if 70^in is unavailable) to perform a x^ fit 
to the Spitzer data. We selected the model which returned the lowest value of x^- The flux 
density at the observed frequency in the best-fitting model SED was then compared with 
the same frequency in the z = model to derive the /c-correction for that target galaxy. 



4. Far and Mid-IR/Radio correlation as a Function of Redshift 

In Fig. la we show the monochromatic 20cm- versus- 70 yum luminosity correlation^. The 
plot shows that galaxies span four orders of magnitude in IR luminosity, and that a small 
number of radio-dominant AGNs are seen above the main trend. However, because of the 
distance-squared stretching of points in this diagram, we prefer to discuss the correlation and 
its dependence on redshift in term of the more basic observable parameters-namely Qir, where 
qir = log (S'ir/S 20cm.), and Sir is the flux density in the 24^m or 70;um bands respectively. 
These monochromatic values for q represent the slope of the IR/radio correlation, and the 
dispersion about the mean q is indicative of the strength of the correlation.^ 

Fig. lb and c show the distribution of qm points before and after the application of a 
/c-correction using the SED model-fitting method. Ignoring the points below q = 1.6 which 
are associated with radio-loud AGN, the figure shows a tight correlation between radio and 
70 /xm emission over a wide range of redshifts extending to z > 1. 

We estimated the median and dispersion of the correlation using a biweight-estimator 
discussed by Beers, Flynn & Gebhardt (1990). This method is resistant to outliers and robust 



"^Here the emitted luminosity is L^^— 47r d;^ S^ofcs/(l+z) (Hogg et al. 2002), where Si,obs has aheady 
been k-corrected for an assumed SED as previously discussed. The luminosity distance, d;, assumes Hq = 70 
km s-i Mpc-i, flM = 0.3, r^A = 0.7. 

^Ideally, it would be better to measure a bolometric q, but insufficient data is available at longer wave- 
lengths to do this reliably (see Papovich & Bell 2002). In addition, monochromatic q relationships may be 
of great interest for studies involving only one or two MIPS bands. 
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Fig. 2. — ^24 as a function of redshift for a) uncorrected and b) k-corrected using the SED- 
fitting method, and c) correcting with an M82 template (see text). 

for a broad range of non-Gaussian underlying populations. The central location (mean) value 
and scale (dispersion) in Qiq before (2.16±0.17) and after (2.15±0.16) correction appears to 
be invariant with 2; up to 1. The slight increase in q^Q above z = 1 is probably the result of 
small number statistics, and the strong Malmquist bias which leads to very luminous galaxies 
appearing in the sample at these redshifts. The highest redshift galaxies have luminosities 
of comparible with that of ULRGS. 

We show in Fig. 2a and b the corresponding q diagrams for the 24/im band. The 
filled circles indicate galaxies that have a measured flux at 70/im. The dispersion in the 
(724- values at zero redshift (mean = 0.84±0.28) is a factor of 1.6 x larger than that seen 
at 70yum. A general flaring of the 524 distribution with z is noticeable in the uncorrected 
data. Fig. 2b, which has been corrected using the SED-fitting method, has fewer points 
than Fig. 2a because the method relies on having detections at 70yum. Since the flaring is 
reduced after the correction, at least part of the higher-z dispersion can be attributed to 
variations in SED, as reflected in the statistics (mean = 0.94±0.23). The deviations from a 
constant q for the highest-z galaxies may be due to Malquist bias which would favor very 
luminousity ULIRGs which may show free-free absorption at radio wavelengths. Fig. 2c 
contrasts the more definitive SED-fitting method with the M82-correction method applied 
to the complete set of 24/im detected galaxies. The formal scatter is also shghtly reduced 
(mean = 1.00±0.27), especially in the core of the distribution, although outlying points are 
less well corrected, presumably because M82 is not a good fit to all the SEDs (formally the 
M82 method reduces the scatter as well as the SED method if the smaller number of points 
only are considered, with the mean = 1.07±0.21). 

The poor 24/im correlations at low z < 0.2 must reflect intrinsic variations in the 
radio/IR ratio within the population. Insight into this can be seen in Fig. 3, which shows ^24 
as a function of the 24/70/im color. This shows that, for a given color, the galaxies have a 
much smaller q'-dispersion. For example, it is clear that the majority of the outlying points 
in the g-distribution are at the extremes of the color distribution, either having very low 
(<0.04 S24/S70) or very high (> 0.2) flux ratios. Note that the AGNs do not follow the clear 
trend with color, allowing easy separation of this population from the star forming systems. 

Also plotted on Fig. 3 are the loci of model a-parameter from the SED library for three 
different redshift ranges. The solid lines represent the model predictions for the full range of 
possible a-parameters represented in the model. The figure shows that the general slope and 
behavior of the g-versus-color diagram is extremely well represented, but the range of possible 
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Fig. 3. — The trend of 0^24 with mid-far IR color. Open circles = galaxies with z < 0.5, and 
filled circles show galaxies with z > + 0.5. The solid lines show the range of q-parameter 
and colors predicted by the DH model SEDs for increasing values of a-parameter. The range 
is from a = 0.063 to 4 (ending with arrow). Curves for three redshift ranges are represented. 

color-space is too limited at present. It will be noted that many galaxies lie just outside the 
range of coverage of the models at low-z, although at high-z the models encompass more 
of the distribution. The consequence of the incomplete coverage of the model parameters 
provides an explanation for the apparent increase in q with z seen in Fig. 2b. The model- 
results are promising, though presently limited in scope. However, they do suggest that 
great improvements in the method will come by the inclusion of more steeply falling SEDs 
in the template library. 



5. Conclusions 

Observations with the Spitzer Space Telescope show that the slope of the 70/um FIR/radio 
correlation is constant to out to ^ = 1, and that the dispersion about the mean is invariant 
to z = 0.5. Further, using the higher sensitivity gained by observations at 24//m , we show 
that, despite the larger scatter, the Mid-IR/radio correlation has constant slope out to 2; = 
2. In a lambda-dominated universe, this implies that the "conspiracy" of finely tuned factors 
which have created the correlation in nearby galaxies existed in the potentially more primi- 
tive galaxies which we observe more than 8-10 Gyrs in the past. Unlike some dwarf galaxies 
seen locally, which are unable to retain their radio-emitting plasma (and deviate from the 
FIR-radio correlation-Klein et al. 1984, Chi & Wolfendale, 1990), these distant galaxies 
must be able to contain their radio-emitting plasma long-enough to fall on the FIR/radio 
relationship. Future work will concentrate on understanding the systematic effects which 
lead to the increased spread in the strength of the MIR correlation, thus allowing the tech- 
nique to be applied to much deeper surveys where significant overlap with powerful sub-mm 
galaxy population will be possible. 
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